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In this paper we explore the set of linear maps sending the set of quantum Gaussian 
states into itself. These maps are in general not positive, a feature which can be 
exploited as a test to check whether a given quantum state belongs to the convex hull 
of Gaussian states (if one of the considered maps sends it into a non positive operator, 
the above state is certified not to belong to the set). Generalizing a result known to 
be valid under the assumption of complete positivity, we provide a characterization 
of these Gaussian-to-Gaussian (not necessarily positive) superoperators in terms of 
their action on the characteristic function of the inputs. For the special case of 
one-mode mappings we also show that any Gaussian-to-Gaussian superoperator can 
be expressed as a concatenation of a phase-space dilatation, followed by the action 
of a completely positive Gaussian channel, possibly composed with a transposition. 
While a similar decomposition is shown to fail in the multi-mode scenario, we prove 
that it still holds at least under the further hypothesis of homogeneous action on the 
covariance matrix. 


1 



I. INTRODUCTION 


Gaussian Bosonic States (GBSs) play a fundamental role in the study of continuous- 
variable (GV) quantum information processin^i® with applications in quantum cryptogra¬ 
phy, quantum computation and quantum communication where they are known to provide 
optimal ensembles for a large class of quantum communication lines (specifically the phase- 
invariant Gaussian Bosonic maps)®^^. GBSs are characterized by the property of having 
Gaussian Wigner quasi-distribution and describe Gibbs states of Hamiltonians which are 
quadratic in the held operators of the system. Further, in quantum optics they include co¬ 
herent, thermal and squeezed states of light and can be easily created via linear amplihcation 
and loss. 

Directly related to the dehnition of GBSs is the notion of Gaussian transformation^®®, 
i.e. superoperators mapping the set 0 of GBSs into itself. In the last two decades, a 
great deal of attention has been devoted to characterizing these objects. In particular the 
community focused on Gaussian Bosonic Ghannels (GBGs)l^, i.e. Gaussian transformations 
which are completely positive (GP) and provide hence the proper mathematical represen¬ 
tation of data-processing and quantum communication procedures which are physically im- 
plemetabld®. On the contrary, less attention has been devoted to the study of Gaussian 
superoperators which are not GP or even non-positive. A typical example of such mappings 
is provided by the phase-space dilatation, which, given the Wigner quasi-distribution Wpir) 
of a state p of u Bosonic modes, yields the function = IF’p(r/A)/A^”' as an output, 

with the real parameter A satisfying the condition |A| > 1. On one hand, when acting on 0 
the mapping 

Wp(r)^IUf)(r), (1) 

always outputs proper (Gaussian) states. Specihcally, given p G 0 one can identify another 
Gaussian density operator p' which admits the function W^\r) as Wigner distribution, i.e. 
lUp(r) = lUj^^(r). On the other hand, there exist inputs p for which IFj^^(r) is no longer 
interpretable as the Wigner quasi-distribution of any quantum state: in this case in fact 
hFj^^(r) is the Wigner quasi-distribution IFg(r) of an operator 6 which is not positiv^® (for 
example, any pure non-Gaussian state has this property for any A 7^ Accordingly 

phase-space dilatations ([^ should be considered as “unphysical” transformations, i.e. map¬ 
pings which do not admit implementations in the laboratory. Still dilatations and similar 
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exotic Gaussian-to-Gaussian mappings turn out to be useful mathematical tools that can 
be employed to characterize the set of states of GV systems in a way which is not dissimilar 
to what happens for positive (but not completely positive) transformations in the analysis 
of entanglement!^. In particular Brocker and Wernei^ used ([^ to study the convex hull € 
of Gaussian states (i.e. the set of density operators p which can be expressed as a convex 
combination of elements of 0). The rationale of this analysis is that the set ^ of density 
operators which are mapped into proper output states by this transformation includes 0 as 
a proper subset, see Fig. Accordingly if a certain input p yields a W^^\r) which is not 
the Wigner distribution of a state, we can conclude that p is not an element of 0. Finding 
mathematical and experimental criteria which help in identifying the boundaries of 0 is 
indeed a timely and important issue which is ultimately related with the characterization of 
non-classical behavior in GV systems, see e.g. Refs.^^^®^, and alscPSEZlfor the fermionic case. 


In this context a classihcation of non-positive Gaussian-to-Gaussian operations is manda¬ 
tory. This analysis has been initiated inpSl^ where Gaussian-to-Gaussian maps are charac¬ 
terized through their Ghoi-Jamiolkowski state, under the hypothesis that this state has a 
Gaussian characteristic function. One goal of this paper is proving this hypothesis: we prove 
that the action of such transformations on the covariance matrix and on the hrst moment 
must be linear, and we write explicitly the transformation properties of the characteristic 


function (Theorem III.l). In the classical case, any probability measure can be written as a 
convex superposition of Dirac deltas, so the convex hull of the Gaussian measures coincides 
with the whole set of measures. A simple consequence of this fact is that a linear trans¬ 
formation sending Gaussian measures into Gaussian (and then positive) measures is always 
positive. Nothing of this holds in the more interesting quantum case, so we focus on it. 


and use Theorem III.l to get a decomposition which, for single-mode operations, shows that 
any linear quantum Gaussian-to-Gaussian transformation can always be decomposed as a 
proper combination of a dilatation ([^ followed by a GP Gaussian mapping plus possibly 
a transposition. We also show that our decomposition theorem applies to the multi-mode 
case, as long as we restrict the analysis to Gaussian transformations which are homogeneous 
at the level of covariance matrix. For completeness we hnally discuss the case of contrac¬ 
tions: these are mappings of the form ([^ with |A| < 1. They are not proper Gaussian 
transformations because they map some Gaussian states into non-positive operators. Still 
some of the results which apply to the dilatations can be extended to this set. 
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FIG. 1. (Color online) Pictorial representation of the structure of the set of states 6 of a CV system. 


^ is the subset of density operators p which have non-negative Wigner distribution (10). 5^ is set 


of states which instead are mapped into proper density operators by an arbitrary dilatation 0 - 
© is the set of Gaussian states and © its convex hull. 6, gl, 5^, and © are closed under convex 
convolution, © is not. For a detailed study of the relations among these sets see Ref.lE!. 


The manuscript is organized as follows. We start in Sec. [TTjby introducing the notation and 
recalling some basic facts about characteristic functions, Wigner distributions and Gaussian 


states. In Sec. HI we state the problem and prove the Theorem III.l characterizing the 
action of Gaussian-to-Gaussian superoperators on the characteristic functions of quantum 
states and its variations, including the probabilistic analog. In subsection III A we consider 


the case of contractions. In Sec. IV we present the main result of the manuscript, i.e. the 
decomposition theorem for single-mode Gaussian-to-Gaussian transformations. The multi- 
mode case is then analyzed in Sec. [Vj The paper ends hence with Sec. ivg where we present 
a brief summary and discuss some possible future developments. In Appendix we prove 
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the unboundedness of phase-space dilatations with respect to the trace norm. 

II. NOTATION 

In this section we introduce the notation and some basic dehnitions which are useful to 
set the problem. 

A. Symplectic structure 

Consider a CV system constituted by n independent modeJ^with quadratures satisfying 
the canonical commutation relations 


Qt, ^ j, A; = . 


( 2 ) 


Organizing these operators in the column vector R of elements 



( 3 ) 


equation ([^ can be cast in the equivalent form 


R\ 1 , 


( 4 ) 


where A is the skew-symmetric matrix 



( 5 ) 


which dehnes the symplectic form of the problem. Accordingly a real matrix S is said to be 
symplectic if it preserves A, i.e. if the following identity holds 


^A^^ = A . 


( 6 ) 


B. Characteristic and Wigner functions 


For any density operator p E & we dehne its symmetrically ordered characteristic function 






This formula makes sense for any trace-class operator p, the associated x(k) being a bounded 
continuous function. By using the Parceval-type formulaP^ 

Tr {p\ P 2 ) = j Xpi(k) Xp2(k) , (8) 

one can extend the correspondence Q to Hilbert-Schmidt operators, the associated x(k) 
being a square-integrable function. The correspondence is the isomorphism between the 
Hilbert space io of Hilbert-Schmidt operators with the inner product given by the left hand 
side of (|^, and that of square-integrable functions of k. The operator p can be expressed as 

P = / xp(k) ^ . (9) 

We also define its Wigner function as the Fourier transform of the characteristic function, 
square integrable for any p E S): 

WM = j x,(k) e--"' . (10) 


C. Heisenberg uncertainty 


The covariance matrix of p is defined as 


with ( R) = Tr 


^ Tr (p ) j | j , (11) 

p R , provided the second moments of p are hnite. The positivity of p 
together with the commutation relations (|^ imply the Robertson-Heisenberg uncertainty 
relation 

a > ±iA , (12) 

where the inequality is meant to hold for both plus and minus signs in the right-hand- 
side. We also recall that the Williamson theoremP^ ensures that given a symmetric positive 
definite matrix a, there exists a symplectic matrix S such that 


SaS^ = ^ua2, 


Z/,; > 0 . 


(13) 


2 = 1 


The z/j are called symplectic eigenvalues of a and can be computed as the positive ordinary 
eigenvalues of the matrix iaA~^ (they come in couples of opposite sign). The Heisenberg 


uncertainty principle (12) is hence satished iff they are all greater than one: 

j = l,...,n . 


Uj > 1 


(14) 
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The symplectic condition ([^ simplifies in the case of one mode. Indeed for any 2 x 2-matrix 
M we have 

MAM^ = AdetM, (15) 

therefore a 2 x 2-matrix S is symplectic iff 

det S' = 1 . (16) 


We recall also that a symmetric positive definite 2 x 2-matrix a has a single symplectic 
eigenvalue, given by 

= det a ; (17) 


then from (14) a > 0 is the covariance matrix of a quantum state iff 


det a > 1 . 


(18) 


D. Gaussian states 


States with positive Wigner function (10) form a convex subset ^ in the space of the 


density operators 6 of the system. The set 0 of Gaussian states is a proper subset of 
A Gaussian state Px,cr ^ 0 with covariance matrix a and first moments x is defined by the 
property of having Gaussian characteristic function 


_ — Ik^crk+ik^x 


x(k) = e J 


(19) 


i.e. Gaussian Wigner function 

W{r) = ^ g-(r-x)^i.-hr-x) _ ^20) 

y det (tt a) 

For a = l 2 n we obtain the family of coherent states Px, ^ ^ IR^"". The Husimi function 
Tr(px, i 2 nP) of bounded operator p uniquely defines p . It follows that the linear span 
of the set of coherent states, and hence of all Gaussian states, is dense in the Hilbert space 
of Hilbert-Schmidt operators i5- Similarly, these linear spans are dense in the Banach space 
of trace-class operators T. 
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E. The Wigner-positive states and the convex hull of Gaussian states 


Starting from the vacuum, devices as simple as beam-splitters combined with one-mode 
squeezers permit (at least in principle) to realize all the elements of 0. Then, choosing 
randomly according to a certain probability distribution which Gaussian state to produce, 
it is in principle possible to realize all the states in the convex hull 0 of the Gaussian ones, 
i.e. all the states p that can be written as 


P= / Px,a dp,{^,a) , 


where p is the associated probability measure of the process. 


( 21 ) 


It is easy to verify that 0 is strictly larger than 0, i.e. there exist states (21) which are 


not Gaussian. On the other hand, one can observe that (21) implies 

1 


IT^(r) = 


3-(r-x) O' hr-x) 


(i/i(x, a) > 0 , 


( 22 ) 


y/det (ttct) 

so also 0 is included into see Fig. [^ . There are however elements of ^ which are not 


contained in 0: for example, any hnite mixture of Fock states 


N 


N 


p = '^p„\n){n\ N<co p„ > 0 ^p„=l 


(23) 


n=0 


n=0 


is not even contained in the weak closure of 0, even if some of them have positive Wigner 
functionl^. 

III. CHARACTERIZATION OF GAUSSIAN-TO-GAUSSIAN MAPS 


Determining whether a given state p belongs to the convex hull 0 of the Gaussian set is 
a difficult problempHSl Then, there comes the need to hnd criteria to certify that p cannot 


be written in the form (21). A possible idea is to consider a non-positive superoperator <F 


sending any Gaussian state into a statd^. By linearity <h will also send any state of 0 into 
a state, therefore if *F(p) is not a state, p cannot be an element of 0: in other words, the 
transformation $ acts as a mathematical probe for 0. In what follows we will focus on those 
probes which are also Gaussian transformations, i.e. which not only send 0 into states, 
but which ensure that the output states $(p) are again elements of 0. Then the following 
characterization theorem holds 







Theorem III.l. Let ^ be a linear hounded map of the space Sj of Hilbert-Schmidt operators, 
sending the set of Gaussian states 0 into itself. Then its action in terms of the characteristic 
function 0 . the first moments and the covariance matrix 0 is of the form 

$ : x(k) ^ X (^^k) , (24) 

$ : X —Kx + yo (25) 

$ : (T —j- KaK^ + a , (26) 

where yo is an vector, and K and a are 2n x 2n real matrices such that a is symmetric, 
and for any a > ±iA 

KaK^ + a> ±iA , (27) 

where the inequalities are meant to hold for both plus and minus signs in the right-hand-sides. 


The condition (27) imposes that *h(p) is a Gaussian state for any Gaussian p. It is weaker 


than the condition which guarantees complete positivity (see Eq. (45) below), which also 
ensures the mapping of Gaussian states into Gaussian states. An example of not completely 
positive mapping fulhlling (27) is provided by the dilatations dehned in Eq. ([^. Such 
mappings in fact, while explicitly not GF^, correspond to the transformations ( [ 2 ^ where 
we set yo = 0 and take 

K = Al 2 „ , a = 0 , (28) 


with |A| > 1 . At the level of the covariance matrices (26), this implies a' = X^a which 


clearly still preserve the Heisenberg inequality (12) (indeed A^a > cr > ±fA), ensuring 
hence the condition ( [2^ . Dilatations are not bounded with respect to the trace norm (see 
Theorem |A.l| of Appendix]^. This explains why Theorem III.l is formulated on the space of 
Hilbert-Schmidt operators. Indeed, via the Parceval formula we can prove that dilatations 
are bounded in this space: 

dk 


ll®(p)ir = J IXp(Ak)r 
= f ixp(k)r 


(27r)’^ 

dk 

(27rA^)^ 


1 

X2n 


For X = X with \p\ > 1 the transformation (28) yields a contraction of the output Wigner 
quasi-distribution. In the Hilbert space the contraction by A is A^” times the adjoint of 
the dilatation by /i = A as follows from the Parceval formula (^. As different from the 
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dilatations, these mappings no longer ensnre that all Ganssian states will be transformed 
into proper density operators. For instance, the vacnnm state is mapped into a non-positive 
operator (this shows in particnlar that the contractions and hence the adjoint dilatations 
are non-positive maps). 


Another example of transformation not fnlhlling the CP reqnirement (45) bnt respect¬ 


ing (27) is the (complete) transposition 


K = T, 


2n 


a = 0 , 


(29) 


that is well-known not to be CP. Unfortnnately, being positive it cannot be nsed to certify 
that a given state is not contained in the convex hnll € of the Ganssian ones. Is there 


anything else? We will prove that for one mode, any channel satisfying (27) can be written 


as a dilatation composed with a completely positive channel, possibly composed with the 
transposition (29), see Fig. We will also show that in the mnlti-mode case this simple 
decomposition does not hold in general; however, it still holds if we restrict to the channels 
that do not add noise, i.e. with a = 0. 

Proof. Let the Ganssian state px.o- be sent into the Ganssian state py^r with covariance 
matrix r(x, a) and hrst moment y(x, a), with the characteristic fnnction 


X-i-(px,A(k) = Xy.r(k) = e 3 


_ ^-Ik^rk+ik^y 


(30) 


We hrst remark that the fnnctions r(x, a) and y(x, a) are continnons. The map <I> is 
bonnded and hence continnons in the Hilbert-Schmidt norm. The reqnired continnity follows 
from 

Lemma 1. The bijection (x, a) —)■ px.o- is bicontinuous in the Hilbert-Schmidt norm. 

The proof of the lemma follows from the Parceval formnla by direct compntation of the 
Ganssian integral 


X/3x,.(k) ,,,(k) 


2 dk 


( 2 vr) 


Next, we have the identity 


Px'jCT'P- x, fT(dx ) Px,a'+a ; 


(31) 


where px, a is Ganssian probability measnre with the hrst moments x and covariance matrix 
a, which is verihed by comparing the qnantnm characteristic fnnctions of both sides. 
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Applying to both sides of this identity the continuous map $ we obtain 


Py(x', cr'), r(x', a') hx, a (dx ) Py(x, r(x, a'+a) 


By taking the quantum characteristic functions of both sides, we obtain 

^ Xy(x',(T'),r(x',(T')(k) /ix,(T(dx) 

Xy(x, (t'+(t), r(x, (k) , k GM. 

Now notice that /ix,o- is the fundamental solution of the diffusion equation: 

du = -dida'^^djU , 

where d is the differential with respect to a, i.e. 


(32) 


(33) 


d=J2 


d 


da'^^ 


and 


R -JL 

' dx< ’ 


(34) 


(36) 


with the sum over the repeated indices. Relation (32) means that for any hxed k, the function 


m(x, a) = Xy(x,(T'+cr),T(x,cr'+cr)(k) is thc solution of the Cauchy problem for the equation (32) 


with the initial condition m(x, 0) = Xy(x,cr'),r(x,(T')(k). Since the last function is bounded 
and continuous, the solution of the Cauchy problem is inhnitely differentiable in (x, a) for 
cr > 0. Substituting 

1 


M(^x, a) = exp 


■ -k^ r (x, cr' + cr) k + ik^ y (x, cr' + cr) 


into (33) and differentiating the exponent, we obtain the identity 

-|k'^ drk + ik^dy = 


= ^ (^^k^diTk- ik'^diY^ dcr*^' (^^k'^djTk- ik'^djY^ 

-^k^ (dida^^djr) k+ ^k^dida^^djy . 

We can now compare the two expressions. Since the left hand side contains only terms at 
most quadratic in k, we get 

diT = 0, (36) 
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i.e. r does not depend on x. Then, the right hand side simplihes into 


- [diY da^Wj-y^) k + dida"-^dj-y . (37) 

Comparing again with the left hand side, we get 

dT{a) = diy da"^djy'^ (38) 

ciy(x, a) = ^dida^^djy . (39) 

Since dr^a) does not depend on x, also dty cannot, i.e. y is a linear function of x: 

y(x, a) = it'(a)x + yo(cr) , (40) 

where K{a) and yo(c’') are still arbitrary functions. But now (39) becomes 

dy(x, a) = 0 , (41) 


i.e. y does not depend on a, i.e. 


y = Tfx + yo , 


with K and yo constant. Finally, (38) becomes 


dT{a) = K da , 


(42) 


(43) 


that can be integrated into 


T{a) = K a + a . (44) 

Thus we get that the transformation rules for the hrst and second moments are given 


by Eqs. (25) and (26). The positivity condition for quantum Gaussian states implies (27) 


The map dehned by (24) correctly reproduces (25) and (26), so it coincides with <F on the 


Gaussian states. Since it is linear and continuous, and the linear span of of Gaussian states 
is dense in it coincides with $ on the whole S^. □ 

Remark 1. A similar argument can be used to prove that any linear positive map <F of 
the Banach space T of trace-class operators, leaving the set of Gaussian states globally 
invariant, has the form (24). By Lemma 2.2.1 oP^ any such map is bounded, and the proof 


of Theorem III.l can be repeated, with S) replaced by T. In addition, since the trace of 


operator is continuous on T, the formula (24) implies preservation of trace. However, the 


positivity condition is difficult to express in terms of the map parameters yo,K,a. 
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On the other hand, if <h is completely positive then the necessary and sufficient condition 


IS 

a > ±i(A — Ak) , (45) 

where 

Ak = KAK^ . (46) 

Thus <h is a quantum Gaussian channeP, and the condition Eq. (27) is replaced by the more 


stringent constraint (45). 

For automorphisms of the G*-algebra of the Canonical Commutation Relations a similar 
characterization, based on a different proof using partial ordering of Gaussian states, was 
hrst given 


Remark 2. There is a counterpart of Theorem III.l in probability theory: 

Theorem III. 2. Let <h he an endomorphism (linear bounded transformation) of the Banach 
space AI(M"') of finite signed Borel measures on IR" (equipped with the total variation norm) 
having the Feller property (the dual 4)* leaves invariant the space of hounded continuous 
functions on IR"/ Then, if $ sends the set of Gaussian probability measures into itself, <h 


is a Markov operator whose action in terms of characteristic functions is of the form (24), 


with the condition (27) replaced by a >0. 


Proof. The proof is parallel to the proof of Theorem III.l, with replacement of (31) by the 
corresponding identity for Gaussian probability measures. As a result, we obtain that the 


action of <h in terms of characteristic functions is given by (24) for any measure /i which is 
a linear combination of Gaussian probability measures. For arbitrary measure fi G AI(M”') 
the characteristic function of *h(/i) is 


X4-(M)(k) - 


A k^> 


$(/r)(dx) = / 




/i(dx) , 


where is continuous bounded function by the Feller property. Since the linear 


span of Gaussian probability measures is dense in A1(M"') in the weak topology dehned by 
continuous bounded functions (it suffices to take Dirac’s deltas, i.e, probability measures 


degenerated at the points of M"") , the formula (24) extends to characteristic function of 
arbitrary hnite signed Borel measure on M”. The action of <h on the moments is given by 


(25) and (26). The positivity of the output covariance matrix when the input is a Dirac 
delta implies a > 0. □ 
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A. Contractions 


A contraction by A = ^ behaves properly on the restricted subset of 0 formed 
by the Gaussian states whose covariance matrix admits symplectic eigenvalues larger than 
Indeed all elements of will be mapped into proper Gaussian output states by the 
contraction (and by linearity also the convex hull of 0 ^( 2 ^ will be mapped into proper output 
density operators). We will prove that any transformation with this property can be written 


as a contraction of l//i, followed by a transformation of the kind of Theorem III.l Let us 
hrst notice that: 


Lemma 2. A set {K, a) satisfies (27) for any a with symplectic eigenvalues greater than /i^ 
iff [fiK, a) satisfies 


for any a > ±iA. 

Proof, a has all the symplectic eigenvalues greater than fj? iO. a > i.e. iff a' = a/pf 


is a state. Then (27) is satished for any a > ±i/i^A iff 

+ a > ±iA \/ a' > AiA , 


i.e. iff {pK, a) satishes (27) for any a > ±iA. 


(47) 


□ 


Then we can state the following result: 


Corollary 1. Any transformation associated with {K,a) satisfying (27) for any state in 


0 


(» 


(i.e. for any a > ±ip‘^A) can be written as a contraction of 1/p, followed by a 


transformation satisfying ([^ for any state in 0 (i.e. for any a > ±iA). 


IV. ONE MODE 


Here we will give a complete classihcation of all the one-mode maps (24) satisfying (27). 
We will need the following 


Lemma 3. A set (K, a) satisfies (27) iff 


Vdeta > 1 — I det K\ . 


(48) 


Proof. For one mode, cr > 0 satishes a > ±iA iff deter > 1, and condition (27) can be 
rewritten as 

det (^KaK^ + «)>!, V cr > 0, det cr > 1 . (49) 
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FIG. 2. (Color online) Pictorial representation of the decomposition of a generic (not necessarily 
positive) Gaussian single-mode transformation in terms of a dilatation, CP mapping and (possi¬ 
bly) a transposition. The same decomposition applies also to the case of n-mode transformations 
when no extra noise is added to the system, see Sec. |V} 


To prove (49) (48) consider first the case det ^ 0. Choosing o such that KaK^ = 

J^g=Q!, we have a > 0, det a > 1. Inserting this into (49), we obtain 


1 + 


I det K\ 

Vdet 


a 


or, taking square root, 


1 + 


I det K\ 
\/det 


a 


det a > 1 


Vdet a > 1 . 


hence (48) follows. 

If det = 0, then there is a unit vector e such that Kg, = 0. Choose o = e“^ee^ -|- eeie^, 
where e > 0, and ei is a unit vector orthogonal to e. Then cr > 0, det cr = 1, and KaK^ = eA, 


where A = KeiejK^ > 0. Inserting this into (49), we obtain 

det {eA -|- a) > 1, V e > 0, 


hence (48) follows. 


To prove (48) (49), we use Minkowski’s determinant inequality 


\/ det(74 -|- B) > a/ det A + \/det B V A, 5 > 0 . 


(50) 


We have for all cr > 0, det cr > 1, 

y/det {KaK'^ + a) > 

> |det K\ \/det cr -|- \/det a > 

> |det K\ + det a >1, (51) 
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where in the last step we have used (48). 


□ 


To compare transformations satisfying (48) with CP ones, we need also 


Lemma 4. A set {K, a) characterizes a completely positive transformation (i.e. satisfies 

Vdet a > |1 — det K\ . 


(45); zff 


Proof. For one mode, using (15), 


and (45) becomes 


Ak = KAK^ = det iP A , 


a > ±*(1 — det K)A . 


Recalling (18), for linearity (54) becomes exactly 

det a > (1 — det K'f . 


(52) 

(53) 

(54) 

(55) 

□ 


We recall here that a complete classihcation of single mode CP maps has been provided 

in Refs. 122134] 

We are now ready to prove the main result of this section. 


Theorem IV. 1. Any map <F satisfying (27) can be written as a dilatation possibly composed 


with the transposition, followed by a completely positive map. In more detail, given a pair 


{K, a) satisfying (27), 

al: If 


0 < det K < 1 , 


(56) 


<h is completely positive. 


a2: If 


det K > 1 , (57) 

<h can be written as a phase-space dilatation of parameter A = a/ det K > 1, composed 
with the symplectic transformation given by 

K 


S = 


a/ det K 

composed with the addition of Gaussian noise given by a. 


(58) 
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bl: If 


— 1 < det < 0 , (59) 

$ can be written as a transposition composed with a completely positive map. 


b2: If 


det K < —1 , (60) 

$ can he written as a dilatation of -^/[deF^ composed with the transposition, followed 
by the symplectic transformation given by 

K 


S = 


composed with the addition of Gaussian noise given by a. 


(61) 


Proof, a: Let us start from the case 


al: If 


det K > 0 


0 < det K < 1 , 


(62) 


(63) 


a2: If 


(|48|) and (52) coincide, so d* is completely positive. 

det K > 1 , 


we can write K as 


where 


K = S Vdet KI 2 , 


(64) 

(65) 

,_ ( 66 ) 

V det K 

is symplectic since det S = 1. Then 4) can be written as a dilatation of Vdet K > 
1, followed by the symplectic transformation given by S, composed with the 
addition of the Gaussian noise given by a. 


S = 


K 


b: If 


we can write K as 


det K < 0 , 

K = K'T , 
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(67) 

( 68 ) 











where T is the one-mode transposition 


T = 




and 


det K' = — det K > 0 . 


From (48) we can see that also K' satisfies 


Vdet a > 1 — |det K' 


(69) 


(70) 


(71) 


and we can exploit the classification with positive determinant, ending with the same 
decomposition with the addition of the transposition after (or before, since they com¬ 
mute) the eventual dilatation. 

□ 


V. MULTI-MODE CASE 


In the multi-mode case, a classification as simple as the one of theorem IV. 1 does not 
exist. However, we will prove that if $ does not add any noise, i.e a = 0, the only solution to 


(27) is a dilatation possibly composed with a (total) transposition, followed by a symplectic 


transformation. We will also provide examples that do not fall in any classification like IV. 1 


i.e. that are not composition of a dilatation, possibly followed by a (total) transposition, 
and a completely positive map. 

We will need the following lemma: 


Lemma 5. 

inf wVw = Iw^Awl V w G . (72) 

cr>±iA 

Proof. 

a. Lower bound The lower bound for the LHS is straightforward: for any a > ±iA 
and w G we have 

wVw > iiw^Aw , (73) 


and then 


inf wVw > Iw^Aw 

a->±iA 


(74) 
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b. Upper bound To prove the converse, let 


W = Wi + iW2 , Wj G 


t)2n 


where without lost of generality we assume wi 7^ 0. Then 

wVw = wfawi + w^crw2, |w^Aw| = 2 |wf Aw2| . 

Assume hrst wf Aw2 = e 7^ 0. Then we can introduce the symplectic basis {e^, 
where 


ei = 


Wi 


hi = 


sign(e) W2 


Expressed in this basis the question (72) reduces to the hrst mode, and the inhmum is 
attained by the matrix of the form 


a = 



1 ) 


where an-i is any quantum correlation matrix in the rest n — 1 modes. 

Consider next the case where wfAw2 = 0 and W2 is not proportional to wi. In this 
context we introduce the symplectic basis {e^, ,1, where 


ei = wi , 


62 = W2 


Accordingly the identity (72) reduces to the hrst two modes, and the inhmum is attained 


by the matrices of the form 


e 0 


cr e = 


0 e 


-1 


e 0 
0 e-i 


< 7^-2 , 


where Un -2 is any quantum correlation matrix in the rest n — 2 modes, and e —)■ 0. 

Finally, if W2 = c wi, c G M, we introduce the symplectic basis {e^, where 


ei = wi. The question (72) reduces to the hrst mode, and the inhmum is attained by the 
matrices of the form 



where Un-i is any quantum correlation matrix in the rest n — 1 modes, and e —)■ 0. □ 

A simple consequence of lemma is 
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Lemma 6. Any a satisfying (27) for some K is positive semidefinite. 


Proof. The constraint (27) implies 


cr + k^ak > 0 


(75) 


for any a > ±iA and k G Taking the inf over a > ±iA, and exploiting lemma with 
w = K^k, we get 

k^ak > 0 , (76) 


i.e. a is positive semidehnite. In deriving (76) we have nsed that, since A is antisymmetric, 
kAk^ = 0 for any real k. □ 


The Lemma allows ns to rephrase the problem: indeed, the constraint (27) can be 
written as 

(A'^w)V(iL^w) + w^aw > |w'^Aw| , (77) 

V (T > ±iA, Vw G Taking the inf over a in the LHS we hence get 


|w'*'A^i'w| + w'I'aw > Iw’^Awl , V w G , 


(78) 


with Ak as in Eq. (46). Notice that, as for the complete positivity constraint (45), since K 


enters in (78) only throngh |wlA/^w|, whether given K and a satisfy (27) depends not on 
the entire K bnt only on A^. 


The easiest way to give a general classihcation of the channels satisfying (78) (and then 


(27)) would seem choosing a basis in which A is as in (|^, and then try to put the anti¬ 
symmetric matrix Ak in some canonical form using symplectic transformations preserving 
A. However, the complete classihcation of antisymmetric matrices under symplectic trans¬ 
formations is very involvecP^, and in the multi-mode case the problem simplihes only if we 


consider maps <h that do not add noise, since in this case the constraint (78) rules out almost 


all the equivalence classes. In the general case, we will provide examples showing the other 
possibilities. 

A. No noise 


The main result of this section is the classihcation of the maps <h that do not add noise 


(a = 0) and satisfy (27): 
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Theorem V.l. A map $ with a = 0 satisfying (27) can always be decomposed as a dilata¬ 


tion 1^ possibly composed with the transposition, followed by a symplectic S transforma- 

K = ST nl2n , (79) 


tion: i.e. 


K = S k1 


2n 


or 


with K > 1. 

Proof. With a = 0 and 


W = Wi + IW2 , 


Wj e 


o2n 


(78) becomes 


|wfAi^W2| > |wfAw2| , 


(80) 


(81) 


i.e. all the matrix elements of A^ are in modulus bigger than the corresponding ones of A 
in any basis. In particular, if some matrix element A^ vanishes, also A®-^ must vanish. Let 
us choose a basis in which Ai^ has the canonical form 

1 


Ai^ - 0 


2=1 


-1 


© 02n—r ; 


where 


r = rank A 


K 


For (|M|), in this basis A must be of the form 

Ai 


A = © 


2=1 


-A,; 


© 02n-r , |Ai| < 1 . 


(82) 


(83) 


(84) 


Since A has full rank, there cannot be zeroes in its decomposition, so r must be 2n. 

We will prove that all the eigenvalues A* must be equal. Let us take two eigenvalues A 
and /i, and consider the restriction of A and A^ to the subspace associated to them: 


/ 


A;. = 


V 




/ 


A = 


/ 


A 


-A 


V 






-/i 


(85) 


If we change basis with the rotation matrix 

^ ( cos^ I 2 — sin6*l2 

\ sin^ I2 cos^ I2 
A ^ R/AR^ , Ak ^ R^kR^ , 


( 86 ) 
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Ak remains of the same form, while A acquires off-diagonal elements proportional to A — /i. 


Since for (81) the off-diagonal elements of A must vanish also in the new basis, the only 
possibility is A = /i. Then all the Aj must be equal, and must then be proportional to 
A: 


Ak = —A 


A 


0 < |A| < 1 , 


(87) 


where we have put all the Aj equal to A 7^ 0 (since A is nonsingular they cannot vanish). 

KAK^ = \a , ( 88 ) 

A 


Relation (87) means 


A (v^77)^ = sign(A)A. 

If 0 < A < 1, we can write K as & dilatation of 

1 



composed with a symplectic transformation given by 


(89) 


(90) 


S=\/XK , 


(91) 


i.e. 

K = S Kl2n , SAS^ = A . 

If — 1 < A < 0, since the total transposition T changes the sign of A; 


TAT^ = -A , 


we can write K as a dilatation of 


K = 


1 


composed with T followed by a symplectic transformation; 


K = ST Kl2n , SAS’’ = A , 


(92) 


(93) 


(94) 


(95) 

□ 
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B. Examples with nontrivial decomposition 


If a 7^ 0, a decomposition as simple as the one of theorem V.l does no more exist: here 
we will provide some examples in which the canonical form of is less trivial, and that 
do not fall in any classihcation like the precedent one. Essentially, they are all based on this 
observation: 


Proposition 1. If a is the covariance matrix of a quantum state, i.e. a > AiA, the 


constraint (27) is satisfied by any K. 


Since for one mode the decomposition of theorem IV. 1 holds, we will provide examples 
with two-mode systems. 

We will always consider bases in which 


/ 


A = 


\ 


(96) 


1. Partial transpose 


The hrst example is the partial transpose of the second subsystem, composed with a 
dilatation of x/u and the addition of the covariance matrix of the vacuum as noise: 


K = \/i/ I ^ I , z/ > 0 , a = I4 
T2 


(97) 


In this case we have 


Ak = 


— V 


\ 


-z/ 


/ 


and i{A — Ak) has eigenvalues 


(98) 


± (1 z/) ± (1 - z/) , 


(99) 
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so that one of them is |1 + \ v\ \ > 1, and the complete positivity requirement (45) 


I4 > ±i(A - A 


K 


(100) 


cannot be fulhlled by any i/ 7^ 0. 

We will prove that this map cannot be written as a dilatation, possibly composed with 
the transposition, followed by a completely positive map. Indeed, suppose we can write K 
as 

K = K' AI4 or K = K'Ti\li, A > 1 . (101) 


Then 


^K' = 


is always of the form (98) with 


z/ 




and also the transformation with K' cannot be completely positive. 


( 102 ) 

(103) 


2. Q exchange 

As second example, we take for the added noise a still the covariance matrix of the 
vacuum, and for the matrix K the partial transposition of the hrst mode composed with the 
exchange of and followed by a dilatation of ^/u: 


( 


a = I4 > ±iA , K = \fv 


\ 


, z/ > 0 . 


With this choice. 


/ 


Ak = 


— V 


\ 


-V 


V 


\ 


The transformation is completely positive iff 


(104) 


(105) 


I4 > ±i{A - A 


K 


(106) 
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and since the eigenvalnes of i(A — Ak) are 


± a /1 + z /2 , 


(107) 


(106) is never fnlfilled for any u ^ 0. 


As before, we will prove that this map cannot be written as a dilatation, possibly com¬ 
posed with the transposition, followed by a completely positive map. Indeed, snppose we 
can write K as 

K = K' AI4 or a: = AT' T 4 AI4 , A > 1 . (108) 


Then 


is always of the form (105) with 


^K' = ±p Aii- 




(109) 


( 110 ) 


and also the transformation with K' cannot be completely positive. 


VI. CONCLUSIONS 

In this paper we have explored both at the classical and qnantnm level the set of linear 
transformations sending the set of Ganssian states into itself withont imposing any fnrther 
reqnirement, snch as positivity. We have proved that the action on the covariance matrix 
and on the hrst moment mnst be linear, and we have fonnd the form of the action on the 
characteristic fnnction. Focnsing on the qnantnm case, for one mode we have obtained a 
complete classihcation, stating that the only not CP transformations in the set are actnally 
the total transposition and the dilatations (and their compositions with CP maps). The 
same resnlt holds also in the mnlti-mode scenario, bnt it needs the fnrther hypothesis of 
homogeneons action on the covariance matrix, since we have shown the existence of non- 
homogeneons transformations belonging to the set bnt not falling into onr classihcation. 

Despite the set of qnantnm states that are sent into positive operators by any dilatation 
is known to strictly contain the convex hnll of Ganssian states £ even in the one-mode casd^S^ 
the dilatations are then conhrmed to be (at least in the single mode or in the homogeneons 


action cases) the only transformation in the class (24) that can act as a probe for <t. 
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Appendix A: Unboundedness of dilatations 

Theorem A.l. For any A 7^ ±1 the phase-space dilatation by X is not hounded in the 
Banach space T of trace-class operators. 

Proof. Fix A 7^ ±1, and let 0 be the phase-space dilatation by A. Suppose 0 to be bounded, 

i.e. 


0 A 


<l|0|| 


A 


V A e T 


Let also 


Eq. (Al) implies 


:= (n |0 (|m)(m|) |n) . 


EIL"'’I < lien VmeN 

n=0 


(Al) 


(A 2 ) 


(A 3 ) 


The moment generating function of i^ 


9m{.q) := 


(m) „-inq _ 


1 — T /I — re*'^ 


n=0 


where g G M and 


Dehne 


r : = 


1 — r e"**? V e®'J — r 

A2-1 
A2 + 1 ■ 

1 — r 


• — 


(A 4 ) 


(AS) 


(A 6 ) 


{/mr(l r) 

Let 0 G ^“(IR) be an infinitely differentiable test function with compact support. We must 
then have 

00 

^(j) {am {n - X'^m)) < || 0 ||oo || 0 || . (A 7 ) 


n=0 


Expressed in terms of the Fourier transform of 


0(fc) = / 0(a;) e'^^^ dx , 


(AS) 
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(A7) becomes 


/ poo ^ 

^ 

n=0 


iX'‘mamk —ikamn <" 

27r J - 
< II 0 IIOO ||0|| . 


Since the sum of the integrands is dominated by the integrable function 

l|0|| 


27r 


m 


we can bring the sum inside the integral, getting 


(j){k) gm (a^k) e 


iX^mamk 


Since for any k 


lim {gm{amk) e 


^ < II 0 IIOO ||0|| . 


iX^mamkX _ g^-|- 


(see subsection A1), by the dominated convergence theorem 

A m a-m k 


lim / (j){k) gm (amk) e* 

J-oc 

roo 0 roc 


27r 


e 3 ^ = / (j)(^x) Ai(a;) dx , 
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where Ai(x) is the Airy function. Now we get 

POO 

/ Ai(a:) (j){x) dx < ||0|| ||0|| 


V d. e 


(A9) 


(AlO) 


(All) 


(A12) 


(A13) 


Since the Airy function is continuous and the set of its zeroes has no accumulation points 
(except — 00 ), there exists a sequence of test functions G r G N with H^r-lloo = 1 

approximating sign (Ai(a;)), i.e. such that 


implying 


lim / Ai(a;) (j)r{x) dx= |Ai(a;)| dx = 00 , 


l| 0 || = 00 


(A14) 


(A15) 


□ 
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1. Computation of the limit in (All) 


Here we compute explicitly the limit in (All). It is better to rephrase it in terms of 


q\= Qmk , q 0 


(A16) 


(remember that am ~ Putting together (All), (A4), (A5) and (A 6 ), we have to 

compute 


1 —r /I —re®'^ ,i+zAni+A \ ? 

- 1—T " 


lim . ... 

g-s>o \ 1—re“®'? \ — T 


= 63 . 


(A17) 


The first term on the left-hand-side tends to one. The second term on the left-hand-side 
instead can be treated via Taylor expansion, i.e. 


for g —0. This gives 


g*<? _ 7 - 3(^x — r)^ ^ 


*:3(1-t)3 

,T —re®^ ,■ 1+T ■;3 t(i+t) 

hm I —^- 6 


q^o \ 6 ®^ — T 

/ iq^T{l + T) 

= hm I 1 H- 7 

q^o y 3(1 — r)^ 

i fc3 

= 63 , 


+ O (/) 


fchi-T)3 

g 3 t( 1 +t) 


(A18) 


(A19) 


which proves the identity of (A17). 
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